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ABSTRACT 

H 2 CO ice on dust grains is an important precursor of complex organic molecules (COMs). H 2 CO gas 
can be readily observed in protoplanetary disks and may be used to trace COM chemistry. However, 
its utility as a COM probe is currently limited by a lack of constraints on the relative contributions of 
two different formation pathways: on icy grain-surfaces and in the gas-phase. We use archival ALMA 
observations of the resolved distribution of H 2 CO emission in the disk around the young low-mass star 
DM Tan to assess the relative importance of these formation routes. The observed H 2 CO emission 
has a centrally peaked and radially broad brightness profile (extending out to 500 AU). We compare 
these observations with disk chemistry models with and without grain-surface formation reactions, 
and find that both gas and grain-surface chemistry are necessary to explain the spatial distribution of 
the emission. Gas-phase H 2 CO production is responsible for the observed central peak, while grain- 
surface chemistry is required to reproduce the emission exterior to the CO snowline (where H 2 CO 
mainly forms through the hydrogenation of CO ice before being non-thermally desorbed). These 
observations demonstrate that both gas and grain-surface pathways contribute to the observed H 2 CO 
in disks, and that their relative contributions depend strongly on distance from the host star. 

Subject headings: Protoplanetary Disks, Astrochemistry, Circumstellar Matter, ISM: Molecules, Radio 
Lines: ISM 


1. INTRODUCTION 

The chemical structure of protoplanetary disks is set 
by a combination of gas and grain-surface chemistry, reg¬ 
ulated by the temper ature, density, and radiation prop¬ 
erties of the disk (e.g.|Aikawa et al.| 2002 |). Radiation im¬ 
pinging upon the disk from both the central star and the 
interstellar radiation field results in strong radial and ver¬ 
tical temperature gradients, differentiating the disk into 
three chemical layers: a cold midplane, a hot atmosphere, 
and a warm molecular layer at int ermediate heights. 


(Aikawa et al. 1996 Bergin et al. 2007). Gas-phase chem- 


istry generally dominates above the midplane. Toward 
the cold midplane, and at larger radial distances from the 
star, the sequential freeze-out of volatiles regulates th e 

Zhang et al.|[^013 ). 


chemical structure jQi et al. 2013b 


Freeze-out of molecules depletes gas-phase abundances 
(affecting gas-phase chemistry), and determines where in 


the disk ice chemistry on gra ins can occur (e.g. Garrod 
[eTallpO^IWalsh et al. ]| 201 ^ . 

Gram-surface chemistry is especially important for 
the production of complex organic molecules (GOMs), 
and thus for the prebiotic potential of planetary bod¬ 
ies forming in the disk (e.g. [Garr od et al.||2008| |Garrod| 
|2013[ iMumma fc Gharnleyj 201 l[ ). With one exception 
(GH 3 GN)(Oberg et al.|2015|, GOMs have not been possi¬ 
ble to observe directly in disks. An alternative approach 
to estimating the GOM abundance in disks is to con¬ 
strain the chemistry that produces them, by observing 
smaller organic molecules with a grain-surface formation 
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pathway. 

Of the molecules with predicted grain-surface forma¬ 
tion routes, only H 2 GO has been found to be readily ob- 
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sequential hydrogenation of GO i ce, a reaction which 
is well-studied in the laboratory (iHiraoka et al. [1994 


Watanabe & Kouchi 2002 huchs et al. ' 2009p . Atomic 


hydrogen for these reactions would likely come from ion¬ 
ization of H 2 . Both the freeze-out of GO and the presence 
of H 2 ionization are well established in disks, based on 
observations of GO snowlines and of ions like HGO+ (e.g. 


^71 

:^i 


Qi et al.|201H|2013bl [Mathews et al. 

T997f 


|2013|[Kastner et S. 

A large uncertainty, however, is how efficiently 

' ice is desorbed into the gas-phase at different disk 
locations. 

H 2 GO can also form through gas-phase chemistryj^ 
One important and well studied gas-phase reaction is 
betwe en oxygen atoms and CH 3 (Fockenberg & Preses| 
2002 ), which is suspected to be a major contributor t o 
H 2 GO abundances in dense cores (Guzman et al.|2013). 

A priori it is not clear whether gas or grain-surface re- 
actions or both govern the H 2 GO gas budget in disks, 
and thus the utility of H 2 GO observations as a tracer of 
grain-surface organic chemistry is not well constrained. 
The spatial distribution of H 2 CO should contain infor¬ 
mation on the contributions of these formation mecha¬ 
nisms. As the major gas-phase formation route of H 2 CO 
has an activation barrier, it will contribute more heav¬ 
ily to the H 2 GO abundance at higher densities (higher 
collision frequency) and temperatures, both of which are 
present in the inner disk. Gas-phase formation of H 2 GO 
should therefore result in a centrally peaked abundance 


^ KIDA reaction network ( Wakelam et al.|20I5 l 
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profile. In contrast, grain-surface formation of H 2 CO re¬ 
quires CO to be frozen out on the grain-surfaces, and it 
should ma inly contribute to H 2 CO gas outside the C O 
snow line (jHenning & Semenov 2008 i|Qi et al.||2013a|). 

Based on the qualitative ditterences in abundance pro¬ 
files expected for different formation pathways (in the gas 
or on grain-surfaces), high spatial resolution imaging has 
the potential to constrain the che mical origin of H 2 CO 
in disks. This approach was used by Qi et al. (2013a) us¬ 
ing Sub-Millimeter Array (SMA) observations of H 2 CO 
toward HD 163296. Their observations of an H 2 CO emis¬ 
sion ring outside the CO snowline revealed the presence 
of grain-surface formation. In an earlier Plateau de Bure 
Interferometer (PdBI) study of the protoplanetary disk 
around the T Tauri star DM Tan, a resolved ring-shaped 
distribution of H 2 CO was also interpreted as evidence 
of grain-surface chemistry exterior to the CO snowline 
( Henning fc Semenov||2QQ8 ). 


well-studied (e.g. Guilloteau & Dutrey 1994| Dutrey 

et al. 

2007 

Oberg et al. 2011 

Teague et al. 2015 

) disk 


(sub-)Millimeter Array (ALMA) observations to evaluate 
whether gas-phase chemistry significantly contributes to 
H 2 CO abundances, or whether grain-surface chemistry 
truly dominates. We compare our observations with de¬ 
tailed chemical models of each formation scenario, con¬ 
straining the importance of both gas-phase and grain- 
surface contributions. 


2. OBSERVATIONS 

DM Tan was observed on 2012 October 23 in Band 7 as 
part of the ALMA cycle 0 project 2011.0.00629.S (PI: E. 
Chapillion). The array configuration included 23 anten¬ 
nas with projected baseline lengths between 17.5 and 380 
m (20-430 kA), providing a synthesized beam of 0.65” x 
0.45” (94 X 65 AU at 145 pc) at a PA of ^-29°. The total 
integration time on-source was 35 minutes. The correla¬ 
tor spectral setup consisted of four windows centered at 
339.674, 340.366, 351.527, and 354.264 GHz, each with 
a total bandwidth of 468.75 MHz in 122.07 kHz (0.105 
km s“^) channels. The observed 5 15 - 4 14 transition of 
H 2 CO at 351.769 GHz is located in the third window. 

The quasar J0510+180 was used for phase calibration, 
Gallisto was used for flux calibration, and J0423-013 was 
used as the bandpass calibrator. We self-calibrated the 
data in GAS A version 4.2 using the DM Tan disk con¬ 
tinuum emission. 

We imaged the continuum emission using the multi¬ 
frequency synthesis mode of the GAS A CLEAN task 
from a total bandwidth of 1.875 GHz, minus channels 
with significant line emission. For spectral line imaging, 
we continuum subtracted the visibilities in the uv plane 
and applied a restoring beam of 1.5” x 1.0” (217 x 145 
AU at 145 pc) along the PA of the original beam to ob¬ 
tain a desirable balance of sensitivity and resolution. We 
generated an initial CLEAN mask using a simple model 
of a disk in Keplerian rotation and then modified it to 
more closely fit the observed emission. The achieved rms 
is ^0.16 mJy beam“^ for the continuum and ~10 mJy 
beam“^ per H 2 GO channel. 

3. RESULTS 


3.1. H 2 CO spectral line emission 
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Figure 1. Channel map of H 2 CO 5i5 - 4i4 emission toward 
DM Tan with contours of 8 (Icr) x [3,5,7] mJy beam“^. The 
data were binned to a velocity resolution of 0.31 km s“^ (i.e. 
3 channels). The LSR channel velocity is shown in upper left 
and the synthesized beam in the lower left of each panel. 


Figure [^presents channel maps of the H 2 GO observa¬ 
tions. The emission is consistent with Keplerian rotation 
about the continuum peak at (aJ 2000 = 04^33"^48.7^, 
JJ2000 = 18°10'09.8". The spectrum presented in Figure 
was extracted from the data cube using the CLEAN 
mask and shows a clear double peaked structure, also 
consistent with emission arising from Keplerian rotation. 

Figure displays the dust continuum. We used the 
same restoring beam as applied to the H 2 GO emission 
to enable comparison of dust (representative of the over¬ 
all disk density gradient) and line emissi on profiles. The 


known inner disk dust cav ity of 19 AU (jEspaillat et al. 
|2Q1Q[ Andrews et al.|[2MT] ) is unresolved. An integrated 
moment-U map of the H 2 CO emission is shown in Fig. 
[ 1 > displaying a large H 2 GO emission disk without any 
apparent hole or gap. The emission extends more than 
a full beam-width beyond the expected G O sn owline lo¬ 
cation from chemical models (see Section 3.2). The mo¬ 
ment map was created using the immoments command 
in GAS A, summing over all channels observed to have 
H 2 GO emission (4-8 km s“^). 

To analyze the H 2 GO distribution, we created radial 
profiles (Fig. |^) by deprojecting the observed emission 
along a principal axis of -24.3° with an inclination of 34° 


cipal ax i 

( [Teague et al.|2Q15[) and azimuthally averaging about the 
continuum peak. 'These profiles do not directly translate 
to column density due to optical depth and excitation ef¬ 
fects, but provide a useful visualization. The continuum 
profile is shown in black and the observed H 2 GO emis¬ 
sion in blue. The H 2 GO emission is centrally peaked, but 
exhibits a broad radial distribution compared to the disk 
density profile, with a half-light radius of ^270. 


3.2. Chemical modeling 

To aid interpretation of the observations we ran a small 
set of disk chemistry models. We use the physical disk 
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Figure 2. H 2 CO spectrum, and images of observed dust and H 2 CO emission from the disk around DM Tau. Panel (a): 
Spectrum of H 2 CO emission, integrated over the CLEAN mask for each channel. Panel (h): Continuum map, extracted from 
line-free channels as described in Section]^ Contours denote 0.16 (la) x [10, 40, 160, 640] mJy beam Panel (c): H 2 CO 
emission, integrated from 4-8 km s“^. The dashed red line denotes the expected CO snowline location. Panel (d): Radial 
profiles of continuum and H 2 CO emission, deprojected and azimuthally averaged. The H 2 CO profile is shown in blue with the 
scale to the right, while the continuum is shown in black with the scale to the left. 


3 

< 

N 


0 200 400 600 0 200 400 600 0 200 400 600 

R [AU] 



Figure 3. Model predictions for the H 2 CO gas-phase density structures (top) and column densities (bottom). Models are 
shown with (middle column) and without (left) grain-surface chemistry. We also show predictions for CH 3 OH (right column), 
an almost pure grain-surface chemistry product. Top: Calculated molecular abundances at various (R,Z) positions within the 
disk. Bottom: Integrated column density as a function of radius. 


the time-dep endent chemical model presented in |Cleev^ 


model of DM Tau derived in I Andrews et al.l (DoTTl) and let al. 119901), but since this is not directly constrained for 

Dk 


et al. (2014), assuming 1 Myr of chemical evolution. The 


initial chemical abundances of the model are listed in Ta¬ 
ble 1. The high energy stellar FUV and X-ray radiation 
field within the disk is calcula ted using the Monte Carl o 
cod e and cross se c tions from Bethell fc BerginI (|2011|), 
see 


Cleeves et al. (2013 


2014) tor details. T'he stellar 
UV spectrum i s taken to be the observed TW Hya FUV 
spectrum from Herczeg et al. (2002 2004), scaled to the 
brighter DM Tau hUV luminosity by a factor of 2.16 
(|Yang et al.||2011|). Finally, DM Tau does not have a 
measured X-rav luminosity, so the X-rav upper limit is 

29 


Damiani et al. 


(1995) of Lxr < 4.6 X 10' 


assumed from 

erg s~^ and the quie scent spectral template presented by 


Cleeves et al. (2013). The standard cosmic ray ioni zation 
rate is assumed (e.g. (qr ^ (3 — 7) x 10 


-17 g-1 


Black 


DM Tau, we also run the models with a reduced rate and 
compare the results. 

The central grain-surface reactions in the network in¬ 
clude the formation of H 2 , H 2 O, CH 4 , NH 3 , H 2 CO, 
CH 3 OH and CO 2 . For the grain-surface reactions, we 
assume that 10 % of the time the grain-s urface reaction s 
result in desorption of the product (e.g., Garrod 2013), 
and as this rate has not been directly measured, we also 
compare to the results of a 1% efficiency. In the gas- 
phase, the primary formation route for H 2 CO is O + 
CH 3 , and thus only proceeds in the upper layers of the 
disk where oxygen atoms are present in the gas. These 
oxygen atoms are mainly supplied by the dissociation of 
gas-phase CO via X-ray generated He+. At most loca¬ 
tions in the disk, this particular H 2 CO-producing reac¬ 
tion is more efficient than any other gas-phase formation 
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Table 1 

Chemical model initial abundances 
Species Abundance® Species Abundance® 


H2 

5.00 

X 

10' 

-1 

H20(gr) 

2.50 

X 

10' 

-4 

He 

1.40 

X 

10' 

-1 

N 

2.25 

X 

10' 

-5 

CN 

6.00 

X 

10' 

-8 

H3+ 

1.00 

X 

10' 

-8 

cs 

4.00 

X 

10' 

-9 

SO 

5.00 

X 

10' 

-9 

Si+ 

1.00 

X 

10“ 

11 

s+ 

1.00 

X 

10“ 

11 

Mg+ 

1.00 

X 

10“ 

11 

Fe+ 

1.00 

X 

10“ 

11 

c+ 

1.00 

X 

10' 

-9 

CO 

1.00 

X 

10' 

-4 

N2 

1.00 

X 

10' 

-6 

c 

7.00 

X 

10' 

-7 

NH3 

8.00 

X 

10' 

-8 

HCN 

2.00 

X 

10' 

-8 

HCO+ 

9.00 

X 

10' 

-9 

C2H 

8.00 

X 

10' 

-9 


Abundances are relative to the proton density 


= 2nH2 • 


pathways for H 2 CO by an order of magnitude or more. 
The H 2 CO and CH 3 OH bin ding energies assumed in the 
present models are 2050 K (iGarrod fc Herb's^|2QQ6 ) and 

*2UU7|, respective!}^ ^^ 
rk 


4930 K (Brown & Bolina 
quently, at all heights in t 


Uonse- 
100 AU, 


le disk outside of 
both of these species are only non-thermally desorbed. 
In the surface layers, where zjr > 0.25, direct UV pho- 
todesorption of H 2 CO and CH 3 OH ices are the primary 
non-thermal desorption mechanisms responsible for the 
gas phase abundances. Below this layer, reactive des¬ 
orption of both H 2 CO and CH 3 OH are the dominant re¬ 
turn processes. The full network includes only the stan¬ 
dard isotopes of C, N, O, and H (e.g., no deuterated 
isotopologues) and contains a total of ^ 6000 reactions 
and ^ 600 species. 

For comparison to the observations, we consider three 
different H 2 CO formation scenarios: ( 1 ) gas-phase only; 
(2) both gas-phase and grain-surface; and (3) grain- 
surface reactions only. Calculations for the first two of 
these scenarios were performed by turning on and off the 
grain- surface reactions for H 2 CO in the [Cleeves et al. 
(2013) network. It is not possible to turn off H 2 CO gas- 
phase formation and destruction pathways in the model 
while leaving reactions on the grains because of the close 
coupling between H 2 CO gas-phase chemistry and the 
production of many other chemical species. For this rea¬ 
son, CH 3 O H (which is domina ted by grain-surface for¬ 
mation, e.g. Hidaka et al.||2004 )) was used to simulate a 
pure grain-surface production scenario. Figure shows 
the resultant radial and vertical abundance profiles of 
each model, as well as the integrated column density pro¬ 
files. The gas-phase only model yields a very centrally 
peaked column density profile; the column density drops 
by more than an order of magnitude inside 100 AU. In 
contrast, the grain-surface only model column density is 
almost fiat out to 400 AU followed by a slow decline. 
In the gas and grain-surface formation model, gas-phase 
formation is visible in the inner disk, with a sharp central 
peak, but outside the CO snowline at ^80 AU there is 
a rise in column density due to grain-surface formation. 
As expected, gas-phase formation regulates the H 2 CO 
abundance in the hot and dense inner disk, while the 
inclusion of grain-surface reactions increases abundances 
outside the CO snow line, where CO has condensed out 
onto grains and is available for H 2 CO formation. 


In addition to radial differentiation, the model results 
show that there is vertical structure in the H 2 CO abun¬ 
dances (top panels of Fig. [3). Grain-surface formation of 
H 2 CO mainly enriches the aisk in H 2 CO gas in a banded 
layer where z/r is ^0.25. This can be explained as a 
chemical balancing act, where the grain-surface forma¬ 
tion and desorption of H 2 CO is most efficient at the CO 
“snow-surface” in the presence of X-rays (or other ioniz¬ 
ing radiation). The grain-surface reac tions to form H 2 CO 
and CH 3 OH in Cleeves et al. (2013) rely on atomic hy¬ 
drogen additions, where H formation is almost entirely 
driven by the ionization of H 2 by X-rays from the cool 
central star. Above this snow-surface layer CO is not 
frozen out in high concentrations, and below this layer 
the disk becomes optically thick to X-rays, making the 
grain-surface chemistry inefficient. 

To explore the robustness of these results to the chem¬ 
ical model inputs that have not been directly measured, 
we ran two additional models with a reduced cosmic ray 
ionization rate of (qr ^ ( 1 ) x 10 “^^ s“^ and a reduced 
reactive desorption efficiency of 1 %, respectively. A re¬ 
duced cosmic ray ionization rate, as might be present 
in a disk with a ” T-Taurios phere” generated by stellar 


winds (Cleeves et al. 2013), reduces the total amount 
of grain surface H 2 CO production, particularly in the 
mid-plane where other forms of radiation are not able to 
penetrate. Reactive desorption similarly affects mainly 
midplane abundances, where it is the main ice desorption 
pathway; at higher disk layers UV photodesorption is a 
more important desorption mechanism. Compared to 
the standard chemical model above, the models with re¬ 
duced cosmic ray ionization or reactive desorption result 
in a substantial (order of magnitude) reduction in H 2 CO 
gas in the disk midplane. These abundance reductions in 
the midplane have no significant effect on the emission 
profile, however, as most of the emission comes from the 
relatively unaffected z/r ^0.25 surface rather than the 
midplane. The predicted radial H 2 CO emission profiles 
for the three H 2 CO formation scenarios displayed in Fig 
[^thus seem robust with respect to the assumed cosmic 
ray ionization rate and reactive desorption efficiency. 

3 . 3 . Simulated observations of ehemieal models 

From the modeled molecular abundance structures, we 
calculated emission profiles of H 2 CO 6 15 - 4 i 4 using the 


non-LTE rad iative transfer code LIME (Brinch & Roger 
heijde 2010). A stellar mass of 0.7 (to determine 

line broadening), inclination of 34°, and distance of 145 




pc we re assumed for the radiative transfer (Teague et al. 
20151). The LIME simulations are shown in the top pan- 
Fig. i ALMA observations were then simulated 


in CASA using simobserve^ assuming the antenna con¬ 
figuration and int^ration time of the original observa¬ 
tions (see Section . The original model output over¬ 
predicted the line nuxes by factors of 2-10. To enable 
comparisons of observed and simulated radial profiles, 
the model abundances were uniformly scaled such that 
the peak integrated flux of the model and the observa¬ 
tions match (middle and bottom panels of Eig. |^. The 
model over-predictions are not surprising, since H 2 CO 
and CH 3 OH are at the maximum extent of the complex¬ 
ity considered in the reaction network, and thus there 
are limited reaction pathways away from these species. 
Additionally, as the destruction products of H 2 CO and 
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( |2008| (E„ = 62.45 vs 32.06 and thus the previously 
noted H 2 CO ring may be partially due to excitation ef¬ 
fects. Supporting this hypothesis, LIME simulations of 
our gas and grain model for the 813 - 2 i 2 transition show a 
substantially reduced inner disk component (not shown) 
compared to the 615 - 4 i 4 transition. 

The emission profile we observe in the outer disk of 
DM Tan (i.e. outside the CO snow line) can be explained 
well by grain-surface formation and desorption of H 2 CO. 
However an explanation that accounts for the centrally 
peaked emission line profile also demands a contribu¬ 
tion from gas-phase reactions. The proposed scenario 
where both gas and grain-surface chemistry contribute 
significantly to the disk H 2 CO emission is confirmed by 
chemical modeling. The model results also demonstrate 
that the relative gas and grain-surface chemistry contri¬ 
butions to the observed H 2 CO gas are highly dependent 
on disk radius. 

5. CONCLUSIONS AND FUTURE DIRECTIONS 


Figure 4. Simulations of H2CO emission toward DM Tan, us¬ 
ing the molecular abundance structures from chemical models 
shown in Fig. Top: Non-LTE radiative transfer simulation 
of the 5 i5 - 414 line of H2CO using LIME. Middle: Simu¬ 
lated ALMA observations scaled to the peak observed flux, 
with the synthesized beam shown in the lower left. Bottom: 
Radial emission profiles, where simulated observations have 
been deprojected and azimuthally averaged similarly to the 
data. The simulated emission is shown in blue and observed 
emission shown in black. 


CH 3 OH have a limited chemistry and are not able to 
form more complex species, H 2 CO or CH 3 OH re-creation 
events are likely. Observations have also shown that CO 
abundances in disks may be depleted and time evolving, 
leading to possible ove r-prediction of mo lecules forming 
from CO in the model (Eavre et al.||2013). 

.ine 


The three models exhibit very diiterent emission pro¬ 
files. The gas-phase only model yields compact cen¬ 
trally peaked emission. Both models including grain- 
surface reactions have substantial H 2 CO emission out 
to several hundred AU. When comparing observed and 
simulated radial profiles, it is clear that the gas-phase 
model severely underpredicts emission in the outer disk. 
By contrast, grain-surface formation alone overpredicts 
emission outside 100 AU. The gas and grain-surface 
chemistry model is most similar to the observations, but 
still overpredicts the relative amount of H 2 CO emission 
between 150 and 300 AU compared to the emission from 
the disk center. 


4. DISCUSSION 

The presented ALMA observations show that the 
H 2 CO emission around DM Tan is fairly exten ded, con¬ 
sistent in radial ex tent with the observations of |Henning| 
& Semenov (2008). In contrast to their detected ring 
shape, we hnd that H 2 CO is centrally peaked. The de¬ 
tection of H 2 CO emission from the inner disk i n this 
data set and its absence in Henning & Semenov (2008) 
can be ascribed both to a combination of higher sensitiv- 
ity and better uv coverage in the ALMA data, enabling 
us to recover weaker emission, and the fact that the two 
transitions have different excitation properties. In par¬ 
ticular, the H 2 CO transition analyzed in this work (615 
- 4 i 4 ) has a significantly higher ex citation energy than 
the 3 i 3 - 2 i 2 transition detected by [Henning fc Semenov| 


Based on ALMA observations and detailed chemical 
modeling, we have shown that there are both gas-phase 
and grain-surface contributions to the observed H 2 CO 
gas in DM Tan. This is the first evidence of gas-phase 
formation of H 2 CO toward a protoplanetary disk, and 
implies that H 2 CO cannot be used as a simple tracer of 
grain-surface chemistry without constraining the relative 
contributions of each formation route. 

The exact morphology of the H 2 CO emission will de¬ 
pend sensitively on the ratio of these contributions, but 
will be difficult to recover from observations due to ex- 


citation deg e nerac ies (as demonstrated in the Henning 
& Semenov ( 2QQ8| ) observations) and the currently un¬ 


known vertical abundance structure. Additional ALMA 
observations of multiple H 2 CO transitions could be used 
to provide constraints on excitation conditions, and 
therefore the vertical H 2 CO structure in a disk. Such 
observations would also provide important constraints on 
ice desorption from grain-surfaces and thus how well ice 


and gas abundances are c oupled in disks (e.g. Oberg et al. 
2007 Eayolle et al.|| 2011 ). 


Tmally, we have conclusively shown that some frac¬ 
tion of the H 2 CO around DM Tan must originate from 
hydrogenation of CO on grain-surfaces. We therefore ex¬ 
pect the next hydrogenation produ ct of H 2 CO, methanol 


(CH 3 OH), to be p resent as well (Cuppen et al. 2009 


Walsh et al. 2014). Although it has yet to be detectec 


toward a protoplanetary disk, searching for both C H 3 OH 


and the complex products of its photochemis try (e.g Gar- 
^[2013 Elsila et al. [2007 Thi et aL]|2004 ) clearly 


rod|ZUib Elsila et al. |ZUUI im et al.|ZUU4) clearly war- 
rants further exploration as the capabilities of ALMA 
continue to improve. 
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